Design and Initial Performance of the Askaryan Radio Array 
Prototype EeV Neutrino Detector at the South Pole 
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We report on studies of the viability and sensitivity of the Askaryan Radio Array (ARA), a new initiative to 
develop a Teraton-scale ultra-high energy neutrino detector in deep, radio-transparent ice near Amundsen-Scott 
station at the South Pole. An initial prototype ARA detector system was installed in January 2011, and has 
been operating continuously since then. We report on studies of the background radio noise levels, the radio 
clarity of the ice, and the estimated sensitivity of the planned ARA array given these results, based on the first 
five months of operation. Anthropogenic radio interference in the vicinity of the South Pole currently leads to 
a few-percent loss of data, but no overall effect on the background noise levels, which are dominated by the 
thermal noise floor of the cold polar ice, and galactic noise at lower frequencies. We have also successfully 
detected signals originating from a 2.5 km deep impulse generator at a distance of over 3 km from our prototype 
detector, confirming prior estimates of kilometer- scale attenuation lengths for cold polar ice. These are also the 
first such measurements for propagation over such large slant distances in ice. Based on these data, ARA-37, the 
200 km^ array now under construction, will achieve the highest sensitivity of any planned or existing neutrino 
detector in the 10^^ — 10^^ eV energy range. 



I. INTRODUCTION 

At PeV energies and above (1 PeV = 10^^ eV), ultra- 
high energy (UHE) neutrinos can be most efficiently detected 
in dense, radio-frequency (RF) transparent media via the 
Askaryan effect 1 1 1. The abundant cold ice covering the geo- 
graphic South Pole, with its exceptional RF clarity, has been 
host to several pioneering efforts to develop this approach, in- 
cluding RICE |2| and ANITA |3, 5|, which trace their roots 
back to suggestions made in the mid-1980's O. Within the 
last year, a multi-phased experiment (the Askaryan Radio Ar- 
ray, or 'ARA") has been initiated at the South Pole |7|, de- 
signed to ultimately accumulate hundreds of GZK neutrinos: 
those neutrinos which arise from the cumulative interactions 
of ultra-high energy cosmic rays (UHECRs) throughout the 
universe, specifically through the Greisen-Zatsepin-Kuzmin, 
or GZK process |8, 9|. Combining the drilling and logistics 
experience of IceCube, the in situ RF experience of RICE and 
the sub-nanosecond-scale data acquisition and triggering of 
ANITA, the singular goal of Phase 1 of ARA is to make the 



first definitive observation of the cosmogenic GZK neutrinos 
with a radio antenna array covering a roughly 200 km^ area, 
denoted here as ARA-37. 

In the longer-term, the modularity of the ARA design will 
engender cost-effective growth of the array to the Teraton 
scale necessary to carry out an observatory-class UHE neu- 
trino science program with broad intellectual merit, compris- 
ing basic particle physics, including neutrino flavor and cross- 
section studies in previously unprobed kinematic regimes; as- 
trophysics such as tests of models of galactic evolution, as 
well as cosmic acceleration and propagation models relevant 
to ultra-high-energy cosmic ray (UHECR) observatories such 
as the Pierre Auger and Telescope Array Experiments. 

A critical question for any radio-based instrument deployed 
at the South Pole is whether the local man-made radio inter- 
ference will prevent such an instrument from reaching its sci- 
ence goals. In addition, in situ measurements of the radio 
noise floor, and the radio transmissivity of the ice, as well as 
its propagation characteristics, are necessary to assess the ul- 
timate sensitivity of ARA. 
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To address these issues, we have deployed an initial ARA 
detector station in prototype form, which we also denote the 
ARA Testbed, with much of the functionality planned for the 
ARA-37 design. One of the features of the current ARA de- 
sign is that each element of the large array comprises a stan- 
dalone neutrino detector for its surrounding ice. Although 
multiple- station events are in no way excluded, they are not 
required for either triggering or reconstruction of the interac- 
tion vertex or neutrino direction. This is done to maximize the 
sensitivity of the array at the expense of some angular and en- 
ergy resolution. Once UHE neutrino fluxes and first-order en- 
ergy spectra are established with a discovery instrument, then 
future installations can be further optimized for resolution as 
well as sensitivity. 

In our current design, a single array element, denoted as a 
station here, consists of a cluster with of order 16 embedded 
antennas, deployed up to 200 m deep in several vertical bore- 
holes placed with tens-of-meter horizontal spacing to form a 
small sub-array. We report here on the initial deployment of 
our first station, and results of measurements over the first four 
months of its operation. In all aspects the environment and ice 
quality have exceeded our expectations, and the prospects for 
ARA-37 thus appear excellent. 



II. INSTRUMENT DESIGN 

The ARA stations must be designed to detect a highly 
linearly-polarized RF impulse arriving from any direction, ei- 
ther as an approximate plane wave (for distant sources) or 
a spherical wave for sources that are closer. The individual 
antennas that comprise a station must have response that is 
preferably dipole-like in angular extent, but the antenna ar- 
ray must also be sensitive in two orthogonal polarizations to 
avoid any bias against detecting arbitrary planes of polariza- 
tion, since the RF impulses that arise from the Askaryan effect 
may be observed at an arbitrary angle with respect to the plane 
of polarization. 

Because of practical limitations on drilling, the antennas 
should fit within a 15 cm diameter borehole, and they should 
be at a depth below the firn, the top layer of the ice sheet 
made of compacted snow of density lower than that of solid 
ice, and extending down typically 150 m at the South Pole. 
Antennas within the firn are subject to ray-bending within the 
refractive index gradient which complicates the triggering and 
reconstruction, and also leads to an effective inverted horizon, 
cutting off the detectability of more distance sources in the 
ice. Since South Polar ice has its highest RF clarity in the 
coldest ice, the antenna array response should be optimized 
for detection in the upper 1 to 1.5 km of ice, although deeper 
events will also be seen. Since drilling costs increase dra- 
matically with the required depth of the hole, there is clear 
cost-benefit in making the holes no deeper than needed to 
achieve the science goals. The array must also guard carefully 
against any self-noise radio interference, that is: all data ac- 
quisition, power, and transmission lines must be electromag- 
netically shielded with extreme care. 

In Table I, we list the primary instrument specifications 



chosen for the ARA station design, giving key parameters 
such as the geometry and antenna requirements. The ARA 
prototype instrument shares much of the functionality of the 
planned complete ARA stations, with several notable differ- 
ences, which are enumerated in the Table. 



A. Overall system 

A block diagram of the installed ARA testbed station is 
shown in Fig. [T] A variety of different antennas, both down- 
hole and near-surface, are used. In each case a low-noise 
preamplifier is inserted in close proximity to the antenna to 
minimize insertion loss and associated thermal noise. Highly- 
shielded coaxial cables (Andrew Heliax or Times Microwave 
LMR- series cables) then lead the signals to a secondary re- 
ceiver which boosts the signal with another 40 dB amplifier 
to bring the thermal noise levels (which are several /iV ref- 
erenced to the LNA inputs) to several tens of mV RMS at 
the DAQ system. The signal is split for separate trigger and 
waveform recording paths. For each trigger, the RF wave- 
forms and associated housekeeping (temperatures, threshold 
settings, antenna singles rates, etc.) are recorded locally by 
a single-board computer, then transmitted over twisted pair 
wires via an ethernet modem to a data receiving computer 
at the South Pole station. There the data is repackaged for 
transmission via satellite link to northern hemisphere com- 
puter archives. Future installations will use an optical fiber 
or wireless data link for the data transmission to and from the 
individual stations. 



B. Antennas 

Antenna designs for ARA are highly constrained by the 
need to fit into a borehole in the ice. Currently, ARA has 
demonstrated the ability to successfully drill 15 cm diameter 
boreholes of the required 200 m depth. This hot- water drilling 
operation requires pumping out of the drill water after comple- 
tion of the hole, to avoid freeze-in of ARA equipment, at least 
in the early phases of array construction where recovery of 
downhole instrumentation when necessary helps to improve 
reliability. For our initial season, equipment for pumping the 
holes dry was not available yet, and thus the initial season's 
boreholes were drilled only to a depth of order 30 m typically, 
where pooling of water prevented deeper drilling. However, 
the 15 cm diameter for these holes was maintained. 

An additional issue for borehole antennas where there are 
multiple antennas along the hole, fed by cables from the sur- 
face, is that the conductive coaxial cables must pass through 
upper antennas to reach lower ones. This constraint forces 
the need for a central hole in the antenna to ensure a uni- 
form azimuthal response, and these holes require re-design 
of the antenna feed region for many traditional borehole an- 
tenna designs. With the diameter constraint, antennas for ver- 
tical polarization are straightforward: a fat-dipole or related 
biconical designs will provide very broadband, well-matched 
antennas. With the additional cable-passthrough constraint, it 
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TABLE I: Design specifications for ARA, and the implementation for the prototype station. 



Specified parameter 



ARA 2012++ planned ARA 2010-2011 prototype 



Number of Vpol antennas 

Vpol antenna type 

Vpol antenna bandwidth (MHz) 

Number of Hpol antennas 

Hpol antenna type 

Hpol antenna bandwidth (MHz) 

Number of Surface antennas 

Surface antenna type 

Surface antenna bandwidth (MHz) 

Number of signal boreholes 

Borehole depth (m) 

Vertical antenna configuration 

Vertical antenna pair spacing (m) 

Approximate geometry 

Approximate radius (m) 

Number of calibration antenna boreholes 

Calibration borehole distance from center (m) 

Calibration borehole geometry 

Calibration signal types 

LNA noise figure (K) 

LNA/amplifier dynamic range 

RF amplifier total gain (dB) 



bicone 
150-850 
8 

quad- slotted cylinder 
200-850 
4 

fat dipole 
30-300 
4 

200 

H,V above H,V 
20 

trapezoidal 
10 
3 

40 (2), 750 (1) 
isosceles triangle 
noise and impulse 

< 80 

30:1 

>75 



2 near- surface, 4 in ice 
bicone 
150-850 
2 near- surface, 6 in ice 
bowtie-slotted-cylinder 
250-850 
2 

fat dipole 
30-300 
6 
30 

V or H above H or V 
5 

trapezoidal 
10 
3 
30 

equilateral triangle 
impulse only 
<80 
30:1 
>75 




NOTES: 

A MINI-CIRCUITS ZFSC-3-4-N+ 
A MINI-CIRCUITS ZFSC-2-2-N+ 
A WESTERMODDW-120 



MAIN POWER & DATA 
POWER 

DATA OR RF SIGNALS 
POWER & RF SIGNAL 
RF CALIBRATION SIGNAL 



FIG. 1: Block diagram of the entire ARA prototype system. 



is more difficult to find a design that retains azimuthal sym- 
metry while still passing the signal cable, but we have settled 
on a hollow-center biconical design, where the feed region is 
annular around the passthrough cables, and is fed at multi- 
ple locations along the annulus, with appropriate impedance 
matching. 



For a corresponding antenna with horizontal polarization 
(Hpol) two designs were implemented for testing in the ARA 
testbed: a bowtie-slotted-cylinder (BSC) antenna, and a quad- 
slotted-cylinder (QSC) antenna with internal ferrite loading to 
effectively lower its frequency response. The goal for both 
sets of antennas was to cover a frequency range from about 
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FIG. 2: ARA testbed downhole antennas: left two images, wire-frame bicone Vpol antennas; right two images, bowtie-slotted-cy Under Hpol 
antennas. 



^^^^^ 



Elevation angle Q( 
4> = 0o 




Azimuthal angle (|) 0° 

= 90° 




FIG. 3: Left: Quad-slot cylinder antenna used in one borehole for ARA-testbed. Center: Simulated Gain (dBi) vs. elevation angle ( zero 
degrees is the vertical direction) for three frequencies for the QSC antenna. Right: Simulated Gain (dBi) in the horizontal plane vs. azimuth, 
showing the high degree of uniformity of the QSC azimuthal response. 



150 MHz to 850 MHz. This goal was achieved with the 
Vpol antennas, but the 15 cm diameter borehole constraint 
has proved challenging for the Hpol antennas, both of which 
have difficulty getting frequency response below about 200- 
250 MHz in ice. In addition, the BSC antenna, although it 
was found to have better efficiency than the QSC, suffers from 
some azimuthal asymmetry in its response, and thus the QSC, 
which has uniform azimuthal response, will be used for fu- 
ture ARA stations. In the current testbed station, we have 
primarily used the BSC antennas because of the ease of their 
manufacture for the 201 1 season. Figure|2]shows photographs 
of the wire-frame bicone antennas and the BSCs as they were 
readied for deployment. Fig |3] shows a photo of one of the 
QSC prototypes (only one of the 4 slots is evident), along 
with simulated results for the gain patterns in elevation and 
azimuth, illustrating the uniformity, which was confirmed at 
several angles in laboratory measurements. 

Figures [4] and [5] show the voltage standing wave ratio 
(VSWR), along with the power transmission coefficient for 
the primary borehole antennas used for the ARA-testbed. 



VSWR is related to the complex voltage reflection coefficient 
p of the antenna via the relation 



VSWRiy) = 



|p(v) + l| 
|P(v)-l| 



and the effective power transmission coefficient T (either as a 
receiver or transmitter from antenna duality) is given by 

r(v) = |i-p(v)|2 

and may be thought of as the effective quantum efficiency of 
the antenna vs. frequency V although RF antennas in the VHF 
to UHF range never operate in a photon-noise limited regime. 

In addition to the coupling efficiency of the antennas, the 
other important parameter for RF performance is the antenna 
directivity gain G, often denoted as just gain, and related to 
the effective power collection area of the antenna via the fun- 
damental relation 
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FIG. 4: Measured and predicted Voltage standing wave ratio 
(VSWR, top) and equivalent power transmissivity (bottom) of the 
bicone antenna, with and without the internal ferrite loading on the 
feedline. The ferrite loading is necessary to prevent coupling of RFI 
via stray induced currents the external shield; these results show that 
the presence of the ferrites within the antenna did not degrade the 
performance. 



FIG. 5: Measured and predicted Voltage standing wave ratio 
(VSWR, top) and equivalent power transmissivity (bottom) of the 
bowtie-slotted-cy Under antenna, along with simulations in air and in 
ice. Measurements were also made with a metal bar passing through 
the center of the antenna to confirm its immunity to an internal con- 
ductor. 



for the speed of light c. A^/f is in turn related to the vector 
effective height of the antenna by 

heff = 2^^^ h 

where Rr is the radiation resistance of the antenna, /z is a unit 
vector representing the antenna axis, and r\ = ^Jiijz is the 
impedance of the medium, r| ^ 12071 in free space |[T9l . 
For simple antennas such as dipoles, the induced voltage into 
a matched receiver load is given by V = ^ • heff/2, and for 
our borehole antennas h = z, where z is the local vertical di- 
rection. For our horizontally-polarized antennas which rely 
on coupling to a slot (the Babinet dual to a dipole), the cou- 
pling may either be thought of as a magnetic induction to the 
slot, or as electric coupling to an effective height which is 
toroidal, e.g. h = ^ where ^ is the azimuth unit vector in cylin- 
drical coordinates. In either case, for our broadband antennas 
which have more complex frequency-dependence, detailed es- 
timates of the antenna response wil require an integral over the 
frequency-dependent behavior. For antennas immersed in a 
dielectric, both the speed of light and impedance are affected. 
The dielectric generally lowers the frequency response of the 
antenna compared to air. 

For our biconical antennas, we found the gain to be com- 
mensurate with that expected from a theoretical bicone, with a 
primary mode that behaves as cos^ of the polar angle, and ex- 



tends for about 1 octave in frequency compared to the turn-on 
of the antenna. At higher frequencies, the response becomes 
multi-modal as expected from a dipole used at its higher- 
order resonance regions, but the frequency bandwidth of these 
modes is much larger than for a resonant dipole, due to the 
broadband biconical response. The QSC Hpol antenna gain 
behaves in a very similar manner to the bicones, except that 
its turn-on frequency is higher; this is consistent with the slot 
antennas performing as a Babinet complementary antenna to 
a dipole. 

For the BSC antenna, the response is more like an inverted 
broadband monopole, as shown in Fig. |6] In the top pane, 
the polar angle response is shown. The antenna is oriented 
horizontally for this plot, but in the borehole it is placed verti- 
cally with the high-frequency response peaking down, giving 
the antenna higher sensitivity to upcoming RF at high frequen- 
cies. In azimuth there is also a front- to-back asymmetry which 
leads to better response for incoming RF for the side with the 
open slot. 

Several other antennas were employed in the ARA testbed. 
Several larger discone (Vpol) and batwing (Hpol) designs 
were submerged about 1 meter below the surface, and these 
performed well, but as they require a 40 cm diameter bore- 
hole they are impractical for future use. We also deployed two 
surface fat-dipole (SFD) antennas (Fig. [7] with a frequency 
response extending over a much lower range than our bore- 
hole antennas: 30-300 MHz. These were emplaced to assess 
the environment for detection of geo synchrotron RF emission 
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BOWTIE SLOTTED-CYLINDER, theta vs frequency, phi=0 (theta=0 points to nadir) 

9D 




270 

FREQUENCY, GHZ 
BOWTIE SLOTTED-CYLINDER, phi vs frequencv, theta = 75 (rel. to nadir) 




/SMUTHREL. TOSLOT,deg 

FIG. 6: Top: BSC gain for a range of frequencies, as a function of 
polar angle. The feed point of the antenna is to the left, and the 
antenna is aligned with the x-axis here. Bottom, azimuthal response 
for a similar range of frequencies. 




FIG. 7: Fat dipole used for low-frequency measurements in ARA- 
testbed. The antenna is about 0.5 m in diameter and 4m long, and is 
installed just below the surface in the trench shown. 



impulses have predominantly lower frequency content than 
Askaryan impulses ifTTl , but have RF spectra which can ex- 
tend to hundreds of MHz and thus will be detected in the 
deeper antennas as well. The SFDs in this case can function 
as a way to conclusively tag HAS radio events and distinguish 
them from other RF interference, if the backgrounds are low 
enough in the VHF band. 



C. Preamplifiers and receivers 

The ARA testbed preamplifier module consists of a front- 
end passband and notch filter designed to block any out-of- 
band signals and to notch one strong in-band signal from the 
South Pole communication system at about 450 MHz, and a 
low-noise amplifier (LNA). Any insertion loss in the front- 
end filter causes a loss of signal combined with an increase 
in thermal noise, and this filter was thus carefully designed to 
minimize in-band losses to a fraction of 0.5 dB or less. The 
notch and the out-of-band rejection of the filter are necessary 
to prevent unwanted anthropogenic noise from saturating the 
LNA and driving it into a nonlinear operating range where its 
gain would be compressed. 

The combination of the front-end filter and LNA con- 
tributes about 80K of input-referenced noise to the signal path, 
and the external thermal noise of the ice at an ambient temper- 
ature of 23 OK (the average extends out to a volume of ice of 
order 1 attenuation length in radius), giving a resulting system 
noise of order 31 OK expected. This leads to an input power of 
about -85 dBm for the LNA. Once the signals are amplified by 
the LNA, they are transmitted to the top of the borehole where 
a 2nd stage receiver amplifier with gain of about 40 dB boosts 
the signals again to a range where they can be digitized. 

The gain and noise figure of the preamplifier+receiver sig- 
nal chain was measured for all ARA testbed channels using an 
Agilent HP8970A noise figure meter with an NIST-traceable 
calibrated noise source. These measurements were done with 
the devices under test at an ambient temperature of about - 
15C. This is warmer than the expected temperatures in the ice 
and at the surface during the austral winter, and thus separate 
measurements were made in a freezer at -55C to determine the 
additional temperature scaling for the noise figure and gain. 

The results for the noise figure and gain as a function of fre- 
quency are shown in Fig. [8] The complete system gain is close 
to 80 dB as expected, and the noise figure of the preamplifier 
modules is in the range of 70-90K over most of the band, ex- 
cept for the region of the notch at 450 MHz, and the high edge 
of the band where both cable and low-pass filter attenuation 
begins to dominate. Given that the ambient temperature of 
the ice is still the dominant contribution to the overall system 
noise, these results are well within the specifications. 



D. Trigger, Digitizer, and Data harvesting 



from cosmic ray extensive air showers (EAS) which arrive 
from interactions in the atmosphere above the array. Such RF 



Although it is currently possible to continuously digitize RF 
signals and use these digital waveforms to assess the possible 
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pre- set threshold. To first order the global trigger rate is given 
by: 
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FIG. 8: Total gain and noise figure for ARA prototype preampli- 
fier s+receivers. the two low-frequency antenna LNAs appear to the 
leftmost part of the plot. 



presence of impulsive events of interest, the power require- 
ments for such a streaming trigger are still prohibitive for our 
design goals, which require that future ARA stations be ulti- 
mately able to function at remote locations far removed from 
large-scale power infrastructure. Thus we only digitize wave- 
forms when a separate analog system determines that some 
interesting pattern of RF signals has crossed above an inten- 
sity threshold. 

The ARA-testbed system uses coaxial tunnel diode detec- 
tors to provide a unipolar signal proportional to the instanta- 
neous RF power over a few nanosecond time scale, and these 
pulses are fed into a discriminator implemented via a field- 
programmable gate array (FPGA) using low- voltage differen- 
tial signal comparators. The FPGA then generates digital one- 
shot signals for each of the 16 input antennas, and these digital 
outputs can in turn be assessed via firmware logic to determine 
if a trigger condition is met. For our current system, the pri- 
mary trigger is set so that when any 3 of the 8 primary Vpol 
and Hpol borehole antennas produce signals crossing thresh- 
old within a 100 ns window, the trigger condition is satisfied, 
and all waveforms are digitized and read out by the data acqui- 
sition system (DAQ) which is shown in block diagram form in 

Fig. 121 

The trigger condition of any 3 of the primary 8 borehole 
antennas crossing a pre-set threshold leads to a strong de- 
pendence of the overall global trigger rate on the individual 
threshold-crossing rates, or singles-rates Rsingle-> of each an- 
tenna, which in turn is determined both by the average power 
level in each antenna, and by the value and stability of the 



k:N ■ 



I 

k=3 



Single ^ 



(1) 



and X is the coincidence window, 100 ns in our case, and = 8 
is the number of antennas that contribute to the coincidence. 
From this it is evident that the global rate depends on at least 
the cube of the singles rate. Since we are triggering on fluctu- 
ations in the antenna power which follow an exponential dis- 
tribution, Rsingie ^ ^^v{-^thresh) Ml, and the global trigger 
rate will depend on the cube of an exponential of the voltage 
threshold, or the underlying power level, illustrating the de- 
gree of sensitivity to both thermal noise changes and thresh- 
old stability. For this reason, it is generally advisable to im- 
plement a noise-riding threshold, which servos (at a typical 
cadence below 1 Hz) the voltage thresholds to maintain ap- 
proximately constant Rsingie^ while recording the thresholds 
used. This creates a nearly negligible change in overall sensi- 
tivity but leads to more stable operation. Such a noise-riding 
system was designed into our current trigger, but has not been 
yet implemented in the operating period we describe here. In 
fact we do see evidence of temperature-dependent drifts in 
global trigger rate in our current data, as will be shown in a 
later section. 

Our current DAQ system has a maximum sustained trigger 
rate of 25 Hz due to limitations of the readout system that 
will be improved substantially for future stations. However, 
rates above 10 Hz incur substantial deadtime for the system, 
so we currently set thresholds to keep the rate below a few 
Hz. This trigger level is quite sensitive, however, as shown 



by Fig. 10 Here the trigger efficiency is shown as a function 
of the SNR of the received signal relative to the thermal noise 
level of the system, which constitutes the ultimate noise floor. 
Our current system is 50% efficient for any incoming impulse 
which gives antenna voltage amplitude exceeding 3.5a above 
thermal noise. 

12-bit waveform analog-to-digital conversion is performed 
using the Large Analog Bandwidth Recorder and Digitizer 
with Ordered Readout (LABRADOR), an Application Spe- 
cific Integrated Circuit (ASIC) developed at the Univ. of 
Hawaii |12|. In our implementation of the LABRADOR, 
eight of the channels are sampled at 2 Gigasamples/second 
effectively, with a 1 GHz Nyquist cutoff, and the remaining 
eight channels are sampled at 1 Gsample/sec with a 500 MHz 
Nyquist cutoff. The latter channels are primarily the surface 
and near-surface antennas. Near- future implementations of a 
new ASIC for ARA will extend the high-rate sampling and 
bandwidth to all 16 channels. 

Once digitized, the data is transferred from the 
LABRADOR to a single-board computer and packe- 
tized. Data is then transferred via an ethernet extender over 
1.8 km of twisted pair cable from the testbed DAQ computer 
to a server residing in a South Pole laboratory, where it 
is archived on disk. A selection of data is then packaged 
for transfer via a satellite link to northern hemisphere data 
servers. 
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FIG. 9: Block diagram of the ARA remote data acquisition system. 
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FIG. 10: Overall trigger efficiency vs. the voltage signal-to-noise 
ratio of the input RF impulse. The SNR is measured with respect to 
the RMS receiver voltage for the baseline thermal noise level 



E. In situ Calibration 

To ensure that the triggering and data acquisition system 
are functioning properly, the trigger threshold is set such that 
we trigger at a low but quasi-continuous rate on thermal ac- 



cidentals. This rate is typically 1 Hz or less. In addition, the 
system is software-triggered at 0.5 Hz to provide a continuing 
record of unbiased waveform samples of the RF background. 
Finally, a GPS-synchronized Rubidium clock triggers a local 
calibration pulser that sends ~ 250 ps impulse to antennas in 
a borehole about 30 m radius from the center of the testbed 
array, and this impulse is split and transmitted from both Vpol 
and Hpol antennas to provide an impulse calibration for the 
array. Arrival times for these impulses at each antenna help to 
determine residual inter-antenna delays for the array, and this 
in turn is used to model the local index of refraction of the ice. 



1. Deep IceCube calibration pulser 

To date all results of long-path ice transparency measure- 
ments for the Antarctic ice sheets have been made via trans- 
mission paths that are largely vertical through reflections off 
the basal interface at the bottom of the ice mass |[T6ll . Such 
measurements produce conservative lower limits to the aver- 
age integrated ice attenuation length. If the bottom reflection 
coefficient is known, the resulting attenuation length measure- 
ment can be stated with higher confidence, although it still 
represents an average over a large range of ice temperatures. 
For South Pole ice, experiments in this vein have yielded sev- 
eral results on the ice attenuation that indicate km-scale atten- 
uation lengths in the upper km of ice 1.13. - 16J . 
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For ARA we are interested not just in near- vertical propaga- 
tion, but propagation over a full range of angles, since neutrino 
interactions can take place throughout the ice target mass. To 
begin to address measurements that will support horizontal or 
slant propagation, we were able in the 2010-201 1 season, with 
the help and cooperation of the IceCube project, to deploy a 
set of pulse generators and antennas at several locations near 
the ARA edge of the IceCube detector, at ice depths of up to 
2450 m. These pulser systems used several kV-peak impulses 
transmitted through a vertically-polarized biconical antenna. 
Estimates of the signal strength received in the ARA-testbed 
using conservative values for the ice propagation were used to 
size the choice of transmitter amplitude, but with some risk 
that, if the ice were much more transparent than expected, 
the ARA-testbed might see saturation in the received signals. 
This risk was accepted since stations further away from Ice- 
Cube would still receive signals within acceptable amplitude 
ranges. Results of these pulser tests are presented in a later 
section. 



III. INSTRUMENT PERFORMANCE 

The ARA testbed installation was completed in the third 
week of January 201 1, at a location approximately 1.8 km grid 
East of the IceCube detector, and has been operating nearly 
continuously since that time. The instrument's internal Ru- 
bidium clock, stabilized to a GPS receiver, provides a 1 Hz 
pulse that is used to trigger the calibration pulser, and this 
in turns provides a stable repetitive RF trigger which can be 
used to monitor trigger and amplitude stability. The system 
also triggers on fluctuations of ambient RF thermal noise at 
quiescent rates of typically 0.5-1 Hz, and in addition we use 
the DAQ computer to force a trigger at 0.5 Hz to get unbiased 
waveforms of the ambient RF noise in each antenna channel. 
During the early part of the operation of the testbed in mid- 
January, we also took several dedicated runs where we turned 
on the IceCube RF pulsers; since that time we have not re- 
peated these runs as it requires the presence of personnel who 
are not available during the winter months at South Pole. In 
the remainder of this section we report on various measures of 
the instrument performance for the first 4 months of 201 1. 




FIG. 11: Top: Recorded South Pole wind speeds near the ARA- 
testbed site. Upper middle: Air temperature near the site. Lower 
middle: internal temperature with the data acquisition electronics 
housing. Bottom: Global trigger rate for the ARA-testbed, showing 
a strong correlation with the electronics temperature. 



after it was deployed, to avoid having any remaining conduc- 
tive structure on the surface. We have so far found no evident 
trigger-rate correlation with wind speeds. 

In the current electronics version, the digital-to-analog 
converter (DAC) voltages that set the discriminator thresh- 
olds are held at fixed digital values, but due to temperature- 
dependence, the actual output voltages of the discriminator 
thresholds will drift slightly with temperature. Our opera- 
tional software has the capability to allow the threshold to 
servo on the overall rates, but this was not implemented for 
the testbed yet, and thus we expect some correlation of the 
overall trigger rate with internal temperature, and in fact we 
see such correlation clearly in the data. 



A. Environmental 



B. Thermal noise floor 



Fig. 11 shows a time series of several measured environ- 
mental conditions for the ARA testbed. We monitored the 
wind speed and outside temperature at locations near the 



South Pole station. Fig. 1 1 shows that the overall average trig- 
ger rate of our system has remained quite low for the duration 
of its operation so far, and shows no obvious correlation with 
wind speed. Possible correlation of RF thermal noise with 
blowing snow was a hypothesis we wished to check, since 
blowing snow is known to create static charge buildup on ca- 
pacitive metal structures in polar regions where the humidity 
is extremely low. For the ARA-testbed, all surface equipment 
was deployed so that snow would drift in and bury it soon 



An in situ estimate of the thermal noise floor for each of the 
antennas and receiver systems was made during the installa- 
tion using a bootstrap method. As noted above, the receiver 
noise figures were measured separately, along with antenna 
transmissivity. During installation, the antenna was replaced 
by a 50Q. load termination at the ambient temperature of the 
field site, and a series of noise sample waveforms were col- 
lected for the load as well as with the antenna replaced imme- 
diately afterwards. Because the termination gives only pure 
thermal Johnson noise compared to the antenna, any change 
in the waveforms RMS noise once the antenna is installed 
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FIG. 13: Average noise power spectral density for the surface anten- 
nas with extended low frequency response. 



is due only to a change in the relative noise power, and this 
gives the antenna temperature as compared to the load tem- 
perature. In practice, the antenna temperatures observed in 
this calibration were all comparable to the load temperature 
within the standard errors of the measurement, which were of 
order 10%. This shows that the antenna temperatures are con- 
sistent with the ambient ice temperature. An exception in this 
procedure was observed for the two surface antennas, which 
were subject to the strong Galactic noise component below 
100 MHz. Above 100 MHz the observed noise in these anten- 
nas was consistent with the thermal noise of the ice however. 



Figure [12] shows average Fourier power spectral density of 
the antenna+receiver thermal noise, taken for several hundred 
events which were selected as unbiased noise sample wave- 



forms from data taken in late April 201 1, run 2533. The ther- 
mal noise power level has been first-order calibrated, with a 
systematic error of about ±1 dB. Pure thermal noise at 290 K 
produces a noise power spectral density of -174 dBm/Hz, 
and in our case our average ice+receiver thermal noise is just 
above this, at about 325 K, about 0.5 dB above room temper- 
ature equivalent. In most cases the data match these expecta- 
tions well, but in some cases our gain calibration appears to 
be offset. 

There are several features apparent in these data that are 
worth noting. First, for each channel the passband filter re- 
sponse is evident. For the borehole antennas (upper eight), 
the passband is 130-850 MHz; for the next 6 of 8, the pass- 
band is 100-400 MHz; these antennas are near-surface anten- 
nas that are physically larger and thus respond down to lower 
frequencies, although their beam patterns are still designed to 
primarily view down into the ice. The last row contains two 
borehole Hpol antennas, the QSCs, and the two low-frequency 
surface fat dipole antennas. The QSCs do not turn on until 
about 200 MHz. Below their turn-on frequency they act as 
effective terminations, so the thermal noise power transitions 
over to Johnson noise. For the surface dipoles, their LNAs 
are ineffective above about 300 MHz, so they are lowpassed 
above that frequency. 

Second, some channels show ripple of up to ±1 dB ( 
~ ±10% in power) in the noise spectra; these are likely due 
to residual antenna impedance mismatches, which create low- 
level reflections between the antenna and LNA. Such effects 
can be mitigated with more careful matching which has been 
achieved in other antennas that do not show these effects. In 
practice this spectral ripple leads to very little distortion for 
pulse measurement since the reflection is causally late com- 
pared to the signal. 

Third, in each of the borehole antennas (the upper eight 
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FIG. 14: Modulation of the average RMS voltages observed in the 
two surface low-frequency antennas during one solar day. 



inclination of ~ 63° with respect to the horizon at the South 
Pole, the dipolar response pattern alternately views the hotter 
galactic plane and the cooler polar regions. We have tested for 
this effect for three months of our testbed data, and these are 
folded together according to a sidereal period. This is shown 
in Fig. 



14 for the two surface antennas, and there is clear si- 



nusoidal modulation evident in the total power, which in this 
case has been low-pass filtered below 70 MHz to enhance the 
galactic noise component. The nominal azimuthal orientation 
of these antennas differs by about 22° according to our as-built 
surveys, and we thus expect a phase shift between the two si- 
nusoids. In fact we observe a 36° phase difference, somewhat 
different than our expectations. However, given that there are 
coaxial cables running within several meters of one of the an- 
tennas, well- within the Fresnel zone of wavefront disturbance 
for these frequencies, we attribute the difference to a distor- 
tion of the main beam of one of the antennas. These results 
in any case lend confidence to our estimates of the system and 
antenna temperatures for the ARA-37 testbed. 



spectra), a dashed vertical line marks the location of the notch 
filter, which is weakly visible in these spectra since some of 
the antenna power is blocked at 450 MHz frequency. In the 
lower antenna channels, the receivers had a 400 MHz lowpass 
filter so the spectra are plotted only up to the point where they 
turn off. In the two surface antennas (last two panes at lower 
left) there is another notch at 46.28 MHz to suppress a known 
radar frequency; this is also evident in the data. For these 
antennas, the rapid rise of Galactic noise below 50 MHz is 
also evident. 

a. Galactic Noise The average galactic radio noise be- 
gins to exceed that of the ambient ice at roughly 150 MHz, 
rising quickly with decreasing frequency V as a power law 
with 



7;;^ = 800K(v/100 MHz) 



-2.5 



(2) 



Fig 13 shows the power spectral density for the average of 
the two low-frequency antennas, showing the excess due to 
Galactic noise at frequencies below 100 MHz. At 60 MHz, 
for example, the average sky temperature is about 3000 K, 
and we observe a noise temperature about 6 dB above our 
thermal noise. For these antennas, the LNAs had a receiver 
noise figure of about 100 K in this spectral range, and the 
dipole response of the antenna effectively averages the noise 
power from the sky above with that of the ice below with 
equal weight. Thus we expect an antenna temperature of about 
1700 K, a factor of about 5 above the antenna+system temper- 
ature due to the ice and LNA, or 7 dB, consistent with expecta- 
tions for galactic noise. The observed slope is also consistent 
with galactic noise; at 100 MHz we see an increase of about 
a factor of 2 over ambient ice+system temperature, consistent 
with equation [2] above . 

Because the two surface low-frequency antenna are hori- 
zontal dipoles, we can make use of the standard dipolar cosine 
response of the antenna beam pattern to test for modulation of 
the galactic noise on the time scale of a sidereal day. As the 
dipole beam rotates through the galactic plane, which has an 



1. Solar Radio burst of Feb. 13, 2011 

On February 13, 2011, the Green Bank Solar Radio Burst 
Spectrometer detected a series of strong type II solar radio 
bursts in a broadband measurement spanning 10 MHz up to 
1 GHz fTS], as seen Fig. [15] (Bottom panels). Fig. [BfTop) 
shows the resulting ARA dynamic spectra for several of our 
antennas for the period of the burst. It is evident that the ARA 
testbed saw several components of the burst sequence clearly 
in several different frequency bands. The observed spectra in 
ARA must be deconvolved from the frequency-dependent re- 
sponse of the antennas to produce dynamic spectra with an 
absolute flux density calibration, and this is outside the scope 
of the current report, but the results demonstrate the sensitiv- 
ity of the ARA system, and also the low background noise 
environment of the South Pole ARA site. Further work is un- 
derway to produce interferometric images of the Sun during 
this event. Fortunately, such events are relatively rare and will 
not pose a significant noise background to ARA observations. 



C. Radio-frequency Interference 

As we have noted above, understanding the radio-frequency 
interference environment of the South Pole locale is critical 
to the effective design and operation of a radio array in its 
vicinity. The ARA-testbed was deployed in the latter part of 
the austral summer, and operations at Amundsen-Scott station 
continued at full pace for some weeks after we began taking 
data, giving us a full exposure to the types of interference that 
are common during the austral summer season there, when 
most major science activities are conducted. 



Figure 16 shows results of the trigger rates as a function 
of time for the ARA-testbed over the entire quarter of opera- 
tion. Not surprisingly, the austral summer shows higher levels 
of sporadic interference. There are several specific periods 
where interference is strong enough to significantly impact 
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FIG. 15: Solar radio burst of Feb. 13, 2011 as observed by ARA 
and Green Bank Solar Rado Burst Spectrometer (GBSRBS). Top: 
dynamic spectra of relative RF power in several of the ARA antennas 
as a function of UT for that day (arbitrary units); Bottom panels: 
dynamic spectra and time series plots of the radio intensity seen at 
Green Bank. 



our operation: during meteorological balloon launches, which 
are conducted twice daily until early March, when we observe 
the frequency reduced to once daily. These balloon launches 
utilize a ^ 400 MHz transponder for data telemetry, and this 
produces strong interference for us for about 1/2 hour in each 
flight. During that time, our system event rate is saturated, 
and the deadtime is nearly 100%, but this leads to a total loss 
of livetime of order 5% , and half of this in the later season 
where the launches are reduced to one per day. 

The other significant cause of strong interference is incom- 
ing or departing aircraft which use a 129.3 MHz communi- 
cation channel. In our current system, while our filters are 
rolling off at that frequency, they are not strong enough to 
suppress this signal below the level where it causes a high 
trigger rate, and again we experience 100% deadtime for a pe- 
riod of typically 20-30 minutes around the time of arrival or 
departure of the flights. During peak activity periods, several 
flights a day may land or depart, leading to an additional 5% 
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FIG. 16: Top: Hourly trigger rates for ARA prototype for about 3 
months of operation of 2011. White space in the plot indicates re- 
gions where stored South Pole data could not be transferred yet due 
to satellite telemetry bandwidth restrictions. Bottom: the daily pro- 
file of triggers, showing the peaks associated with weather balloon 
flights. 



per day loss of livetime on average. 

Other than these two sources, there are occasional sporadic 
interference episodes associated with other activity, but these 
have had only a minimal effect on the livetime to date. 



Figure [16] (top) also shows how the trigger rate has slowly 
declined since the summer peak; this is likely not due to de- 
creasing interference, but rather the temperature dependence 
of the trigger thresholds, as discussed above. Overall the 
RFI environment in the locale near Amundsen-Scott station 
is quite acceptable for our requirements, and we estimate the 
yearly average loss of livetime due to interference at of order 
3% at most. 

In addition to assessing noise generated by equipment al- 
ready in place at the station, we installed a wind turbine sys- 
tem at a distance of order 800 m from the ARA-testbed as 
part of a program to investigate autonomous long-term power 
sources (detailed in Appendix A) for ARA-37 and future ar- 
rays that are deployed further out from the existing South Pole 
infrastructure. We have searched our data for triggers which 
specifically reconstruct toward the direction of the wind tur- 
bine, and we find none, indicating that this technology will 
likely be acceptable for use for ARA. 
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D. Calibration Pulser 

The local calibration pulser, which pulses at a repetition 
rate of 1 Hz during all periods when the ARA-testbed sys- 
tem is running, provides a continuous measure of the trigger 
functionality, and a stable timing and amplitude reference for 
monitoring of the ongoing performance of the system. An 
example of an individual calibration pulser event is shown in 
Fig.[T7] Because the calibration pulser is located at nearly the 
same depth as most of the borehole antennas, the signal shape 
as seen by those antennas has the broadest bandwidth and thus 
tends to be very narrow and coherent in time. Because the ra- 
dius of the calibration pulser borehole is only three times the 
~ 10 m radius of the inner borehole array, an amplitude gra- 
dient across the array is quite evident in the data. 

We also noted in these data that there is a second pulse, 
which we hypothesize is a reflection, appearing in the Hpol 
antennas only. Other tests done during deployment using dif- 
ferent calibration locations saw no similar reflections, so it 
does not appear to be something intrinsic to these channels. 
The region where we deployed the ARA-testbed is still rela- 
tively close to the original station established in the 1950's. 
The depth at which we deployed our antennas, typically 20- 
30 m, is in fact the depth at which the 1950 surface layer is 
found. Because many air-drops of equipment in the vicinity 
of the Pole were done during this period, and much equipment 
was left in the field, it is possible we have detected a horizon- 
tal metallic object within several meters of our boreholes, and 
we are investigating this further in the data. 








FIG. 17: Local calibration pulser event waveforms in all ARA- 
testbed antennas. 



One of the important functions of the calibration pulsers is 
to provide a stable reference signal that can be used to de- 
termine the geometry of the array in situ, via the relative de- 
lays between the received waveforms within the array. Fig. [18 
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FIG. 18: Example of delay histogram from cross-correlated wave- 
forms for an ensemble of about 2600 calibration pulser events, in 
this case for borehole antennas VI and V4. The origin is arbitrary in 
this case. 



shows a typical example of a histogram of the delays between 
Vpol borehole antennas VI and V4 for an ensemble of about 
2600 calibration pulser events. The delays for each event 
are extracted by cross-correlating the received waveforms for 
each event, and this delay distribution is shown in the plot, 
with a fitted Gaussian to the peak region giving a well-defined 
mean value (arbitrary in this case) with an estimated error of 
about 4.5 ps, and a Gaussian standard deviation of 136 ps. The 
standard error on the fitted mean represents the ultimate limit 
of the calibration pulser data for determining geometry within 
the array - in this case the 4.5 ps resolution corresponds to 
about 1 mm in ice. Thus the timing resolution is more than 
adequate for calibration; our reconstruction requirements for 
the geometry can easily tolerate few cm accuracy on locations. 

The standard deviation of the fitted Gaussian for the pulser 
delay distribution also provides a calibration of the potential 
angular resolution for the system. A standard deviation of 
136 ps in this case corresponds to about 2.7 cm in the ARA- 
testbed ice. Our angular reconstruction of signal arrival direc- 
tions depends on the ratio of our standard error on the wave- 
form arrival to the length of the projected baseline between 
two antennas, eg: 



AO, 



c X/ 



nSB^ 



ij 



where Xij = ti — tj is the baseline delay for wavefront arrival 
times ti and tj at antennas /, 7, and S is the Poynting vector 
of the incoming wave, here assumed to be a plane wave for 
simplicity. For typical baseline lengths of 15 m or more, the 
implied angular resolution per baseline is of order 0.1° for 
136 ps timing error, again more than adequate for our require- 
ments. 
Fig. 
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shows the distribution of the reconstructed az- 
imuthal direction of events observed on two separate typical 
days during the operation period. Events that RF-triggered the 
testbed during these periods were either thermal noise fluctua- 
tions, which can form random coincidences but do not recon- 
struct as plane- or spherical-waves when analyzed offline, or 
RE triggers from the calibration pulser, or from any other RE 
interference or backgrounds that may be present. The latter 
two types of events recur from fixed directions: the calibra- 
tion pulser from a fixed ~ 30° azimuth; RFI can occur from a 
variety of directions but will be associated with the South Pole 
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FIG. 19: Examples of the azimuth distribution of events observed 
during two days of the ARA-testbed operation. 



tial. Thus we have looked for a background of random triggers 
that are associate with higher-than-normal winds, or structures 
such as our wind-turbine system. 

We find that almost all RF trigger events observed have an 
azimuth consistent with that of the calibration pulser with re- 
spect to the center of the array. January 29, 2011 was a pe- 
riod of relatively high winds (> 20 kt), which produces more 
opportunity for static discharge on any above-surface struc- 
tures, and also yields higher power generation in the nearby 
wind turbine system. By contrast March 19 was a period of 
relatively calm wind. In neither case did we observe any sig- 
nificant number of events that reconstructed to any other di- 
rection, including the ~ 155° azimuth of the nearby wind tur- 
bine system, and other directions associated with the IceCube 
Laboratory and South Pole Station. Fig. [20| shows a two- 



dimensional map of the reconstructed events for the higher 
wind date, with an inset showing the azimuth distribution. A 
very small fraction of events misreconstruct in elevation an- 
gle as evident in the plot, but the vast majority of events are 
consistent with the calibration pulser only. 



E. Deep pulser measurements 

In our measurements, all of the IceCube pulsers show 
some degree of saturation in the received signals at the ARA- 
testbed. While this bodes well for future studies of long- 
range ice transmission for ARA, the saturated waveforms in 
the ARA-testbed complicate efforts at precision estimation ice 
transmission parameters. 




Station, and are likely to repeat, thus giving clusters of events 
in stable directions associated with human activity. It has also 
been hypothesized that wind-blown snow could create elec- 
trical discharges itself, or by depositing charge on structures 
which then discharge once they reach a high-enough poten- 



FIG. 21: Deep calibration pulser event waveforms in all ARA- 
testbed antennas. 



In addition, the bicone antenna design for deployment con- 
centric to the IceCube cables, and under the extreme high 
pressure of the re-freeze that occurs in the hot- water-filled 
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IceCube boreholes, led to some challenges in the antenna con- 
struction that also came with risk. For example, the clam-shell 
jointing of the antennas to surround the IceCube cable must 
be done during a rapid deployment, and inadequate contact 
for the joint can lead to shunt-capacitance across the slot, and 
cause some radiation of Hpol as well as Vpol signals. Also, 
the long-wire vertical IceCube cable in close proximity to the 
antenna can lead to re-radiation and scattering of the signal. 
As a result of analysis of the signals and comparison with 
some prior lab measurements, we conclude that the transmit- 
ted pulse from the deep pulser was dispersed in time over typ- 
ically a ~ 25 ns time window, although it appears that most of 
the pulse power was in fact contained withing this dispersed 
pulse. 

Despite these difficulties, we have observed signals from 
these transmitters at up to a 3.2 km slant range from the 



source. The best example of this is shown in Fig. [21] which 
comes from the furthest and deepest of the deployed pulsers, 
installed in IceCube Hole #1 at a depth of 2450m and a hor- 
izontal distance of about 2 km from the ARA-testbed, for a 
total slant range of 3.16 km at an angle of about 40° from the 
vertical. Some of the observed signals still saturate our DAQ 
dynamic range even at this distance, and we do observe sig- 
nificant Hpol signal as well as Vpol; suggesting that the clam- 
shell coupling of the antenna was not complete. The Hpol 
signals also have higher frequency content than the Vpol sig- 
nals on average, and are factors of 2-3 weaker in amplitude, 
which is consistent with radiation from partial slots along the 
antenna's vertical joint. In several cases, receiving antennas 
are partially shadowed by antennas in the foreground with re- 
spect to the pulse arrival direction. The radius of the first Fres- 
nel zone is the region around any obstruction within which a 
propagating wave is disturbed, and for the ARA-testbed an- 
tennas for a source at distance Rs and the foreground antenna 
at distance d from the receive antenna we have 



Rf 




(3) 



For Rs = 3200 m, and J 20 m, 7?/^ ~ 4 m for ^ Im in 
ice, indicating that we would expect several of the antennas to 
suffer from significant shadowing, which we observe. 



1. Attenuation length estimates 

If we isolate our attention to the foreground Vpol antennas 
that are undisturbed, we may use the resulting amplitude to 
estimate the average attenuation length (La) over the entire 
3.16 km propagation path, along with our knowledge of the 
antenna and pulse parameters at the source. The attenuation 
length can be determined from the Friis equation 1 19], which 
relates the received power to the transmit power Pt as a 
function of frequency: 



Pr{v)=4nPt{v) e 



-2R/{La 



Ar At 



(4) 



where R is the path distance through the ice, and A^,A^ are 
the receive and transmit antenna effective areas (including an- 
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FIG. 22: Top: Antenna individual V2 received waveform from the 
deep pulser, showing the secondary surface reflection near the end of 
the trace. Bottom: cross-covariance of the first pulse with the entire 
waveform. Note that the peak covariance from the second (reflected) 
pulse is inverted in phase with the primary pulse, at a lag time of 
122 ns, which is within 1% of the propagation time estimate for the 
reflected ray. 



tenna coupling efficiency), which are in general also func- 
tions of frequency v. Inverting this for the average attenuation 
length we have: 



-2R (log 



R' 



Pt \Anv^ArAt 



(5) 



The quantity (La) in this case is a weighted average, where 
the weights are determined by the vertical temperature pro- 
file of the ray paths, since the attenuation coefficient (in units 
of m~^) is a = = a(v, Tice), a strong function of the ice 
temperature Tice- 

Accounting for the transmitted power, for the bicone trans- 
mitter antenna response, estimated insertion losses, antenna 
beam pattern, a first-order estimate of the signal saturation, 
and the receive antenna response, our estiamte of the atten- 
uation length over the referenced to 300 MHz is plotted in 



Fig.[23]as a function of depth. The depth dependence is decon- 
volved from the total attenuation under the assumption that the 
depth dependence is dominated by the ice temperature profile, 
which has been determined quite well for South Polar ice. The 
average attenuation length for all depths, 820(+120, —65) m 
is also given. The errors here are dominated by our system- 
atics on the received amplitude and pulse shape. The result- 
ing value is quite consistent with values derived from vertical 
bottom-bounce methods at the South Pole fT6ll . and the sys- 
tematics here are different from those in the bottom-bounce 
measurement in which the bottom reflection coefficient is un- 
known. Note that the attenuation length derived is domi- 
nated by the extended propagation distance in the relatively 
warm basal ice; considering the upper 1.5 km where the ice 
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FIG. 23: The field attenuation length over the 3.16 km path from the 
deep pulser to the ARA testbed as a function of depth. The average 
attenuation length over all depths is also given. 



FIG. 24: Ray paths for direct and reflected rays to antenna V2. The 
red-dash line gives the geometric ray path for the direct ray. 



is much colder than the basal ice, the effective average atten- 
uation length is nearly 1.7 km. This remarkably long field 
attenuation length of the ice in this extremely cold region of 
the cryosphere sets it apart from any other natural (and in 
fact most man-made) dielectrics as the most transparent solid 
medium on Earth. It is perhaps the single most important char- 
acteristic of the ice for enhancing the performance and viabil- 
ity of an array such as ARA. 



2. Surface reflection 

Another interesting feature observed in the deep pulser sig- 
nals, which arrive at the ARA-testbed as a close approxima- 
tion to a plane wave, is the surface reflection which is seen in 
several antennas, most clearly in borehole bicones VI and V2. 
These differ in depths by several meters, leading to a different 
delay time before the reflection arrives. In antenna V2, the 
reflection is clearly phase-inverted compared to the primary 
signal; VI 's record cuts off too early to observe this clearly. 
To verify the presence of the reflection, the waveform (top) 
and a cross-covariance estimate (bottom) of the primary and 
secondary pulses with the primary pulse are shown in Fig. 22 
The cross-covariance has its highest correlation with the pri- 
mary pulse at a 122 ns delay, where the peak is a 60% anti- 
correlation, as seen in the Figure. V2's signal amplitude for 
the reflection is also quite close to the direct signal amplitude, 
which is consistent with a total internal reflection off the top 
of the ice, as expected for the estimated 57° incident angle of 



arrival with respect to the ice surface normal. 



Fig. 24 shows numerical ray trace paths for this antenna 
from the IceCube deep pulser (top), and a zoom to the region 
of the reflection in the bottom pane. The numerical ray tracer 
uses the eikonal equation and a refractive index model devel- 
oped based on radio index of refraction data from the RICE 
experiment at the South Pole, and returns both the ray path 
and total wave transit time along the ray. The delay observed 
in antenna V2's waveform from the direct to inverted reflec- 



tion is 120.0 ns. For the ray paths from Fig. 24 the direct 
path transit time is 18,724.53 ns, and the reflected ray path re- 
quires 18,669.92 ns to the surface, followed by 176.6 ns from 
the surface to the antenna, for a relative delay of 122 ns with 
respect to the direct path, a result identical to that observed 
in the data. This also independently confirms the precision of 
our index-of-refraction model. 

It is a generic feature of a subsurface radio array that signals 
may be observed both directly and via the surface reflection, 
or even via near-reflections for rays that are curved back down 
near the surface; this effect is not yet incorporated into our 
simulations, but will improve the detection efficiency in gen- 
eral, and give some additional reconstruction power as well. 



IV. EXPECTED PERFORMANCE FOR ARA-37 

We characterize the expected performance of ARA based 
on an extensive suite of Monte Carlo and other simulation and 
modeling tools developed over the last two decades. Several 
of the investigators on this proposal have been involved in ra- 
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FIG. 25: Planned layout of the 37 ARA stations with respect to the 
South Pole Station and associated sectors. 
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dio methods for detection of UHE neutrinos since the mid- 
1990's and thus the heritage of our simulation methodology is 
many generations deep and has been proven over a wide range 
of active and completed experiments, including the RICE O 
and AURA [ 20 1 experiments, which are direct pre-cursors to 
ARA, the Goldstone Lunar Ultra-high energy neutrino Exper- 
iment (GLUE) lEB, and the ANITA experiment O. 

Our Monte Carlo tools include detailed ice attenuation 
modeling and raytracing to account for the gradient in the in- 
dex of refraction of the ice vs. depth. They provide state- 
of-the-art modeling of the Askaryan radiation from show- 
ers via tested parametrizations, which have been validated 
by direct measurements of the Askaryan effect in ice at 
SLAC |22|. Neutrino propagation through the earth and ice 
sheets is modeled in detail, and the particle physics of the in- 
teraction, including neutral and charged-current effects, fully- 
mixed neutrino flavors, and secondary shower production due 
to charged-current x- and /i— leptons are accounted for in the 
models. Finally, the detectors are also modeled with high fi- 
delity, including the effects of Rician noise in the detection 
process, spectral response functions of the antennas, and full 
3-D polarization propagation for the radio waves that interact 
with the detector. We thus have reason to report these perfor- 
mance estimates with some confidence. 



A. Station Spacing 

In initial studies of the ARA-37 array |23|, we adopted a 
1.33 km spacing between stations as a compromise which 
gave adequate sensitivity while still allowing for enough over- 
lap between the stations' effective target volumes to yield rea- 
sonable fraction of multi- station coincident events. However, 
in the interim, indications of possible heavier nuclear com- 
position for the UHECRs have led to estimates for the cos- 



FIG. 26: Simulation of the effective acceptance for ARA-37 as a 
function of the inter- station spacing, for an energy of 1 EeV. 



mogenic neutrino flux which are significantly lower than for 
a pure proton spectrum, and thus in our current designs, we 
optimize the ARA-37 array for discovery potential, to maxi- 



mize the number of detected events. Fig. 26 gives the results 
of a parametric study of the neutrino acceptance of the array 
as a function of the spacing for a neutrino energy of 1 EeV, 
which is an excellent proxy on average for the total number of 
detected events integrated over a typical cosmogenic neutrino 
spectrum. 

We find that the acceptance grows rapidly with spacing 
starting with sub-kilometric sizes, and then becomes fully sat- 
urated at about 2-3 km. In our current design we have thus 
adopted 2 km as our baseline spacing; this choice reflects a 
good compromise between maximizing the sensitivity while 
still recognizing the potential logistical costs that grow with 
array size. 



B. Sensitivity 

The primary metric for detection of cosmogenic neutrinos, 
which are presumed to arrive isotropically on the sky, is the 
volumetric acceptance l/Q., in units of km^ steradians. An 
equivalent parametrization is the areal acceptance RQ, (km^ 
sr) and the two are closely related by RQ, = ^Q,/ Lint{Ey) 
where Lint is the interaction length of the neutrinos as a func- 
tion of neutrino energy Ey. The volumetric acceptance, di- 
vided by the instrumented target fiducial volume, gives a mea- 
sure of the detection efficiency of neutrinos which interact 
within the fiducial volume of a detector. In the case of ARA- 
37, one realization of the simulation uses a cylindrical ice 
target volume of radius 10 km, and depth 2 km. Because 
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FIG. 27: Top (large pane): Simulated neutrino volumetric accep- 
tance (km^ sr water equivalent) for the ARA instrument baseline de- 
sign. Bottom left: Depth distribution of simulated events for different 
neutrino energies, showing the contribution of deep ice down to 2 km 
or more at the higher energies. Bottom right: zenith angle distribu- 
tion of detected neutrino arrival directions for a range of neutrino 
energies. Events are detected over a range from ^ 45° above the 
horizon to ~ 5° below it. 



of earth attenuation, neutrinos arrive almost exclusively from 
above the horizon, nominally giving ~ 2% steradians for the 
solid angle. The net target acceptance of the simulated ARA 
detector is thus just over 4000 km^ sr, and this represents 
the maximum neutrino volumetric acceptance the simulation 
could obtain. However, once the constraints given by the neu- 
trino interaction cross section, the ice attenuation length, and 
the ray-tracing geometry of the ice, and the antenna response 
functions are imposed, the effective acceptance becomes an 
energy-dependent fraction of the initial target acceptance. 

shows the simulated VO. results for our adopted 



Fig. 
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baseline design, as a function of neutrino energy in the range 
of interest for the cosmogenic neutrino flux. The acceptance 
reaches the level of > 200 km^ sr at the mid-range of the cos- 
mogenic neutrino flux, which has a broad plateau from about 
5 X 10^^ eV up to just over 10^^ eV, and continues growing 
slowly up to the highest simulated energies, approaching a 
Teraton-steradian. 

Fig. [27jbottom) gives a plot summary of some characteris- 
tics of the simulated data vs. neutrino energy. On the bottom 
left, the depth distribution of detected events is shown nor- 
malized to the event fraction per hundred meters. Events orig- 
inating from below about 2 km depth tend to be suppressed, as 
the attenuation of the ice begins to grow quickly in the warmer 
basal ice l,16J . The estimated average attenuation length of our 



2 km-deep fiducial volume is about 1.5 km, a factor of three 
better than ice in locations such as the Ross Ice Shelf, where 
the thickness is limited to several hundred m, and the attenu- 
ation lengths are comparable to this thickness scale. SP ice, 
especially in the upper 2 km of its depth, is the clearest solid 
dielectric medium on Earth in the radio range, and is the most 
compelling natural feature of the ARA site. 

Fig.[27jbottom) also shows the arrival zenith angular distri- 
bution of neutrino events that were detected, showing that the 
neutrino angular acceptance spans a range from ~ 5° below 
the horizon to ~ 45° above the horizon, more than 6 steradi- 
ans of solid angle. 



TABLE II: Expected numbers of events A^v from several UHE neu- 
trino models, comparing published values from the 2008 ANITA-II 
flight with predicted events for a three-year exposure for ARA-37. 



Model & references Ny-. 


ANITA-II, ARA, 
(2008 flight) 3 years 


Baseline cosmogenic models: 






Protheroe & Johnson 1996 (27) 


0.6 


59 


Engel, Seckel, Stanev 2001 (28) 


0.33 


47 


Kotera,Allard, & Olinto 2010 (29) 


0.5 


59 


Strong source evolution models: 






Engel, Seckel, Stanev 2001 (28) 


1.0 


148 


Kalasheveffl/. 2002 (30) 


5.8 


146 


Barger, Ruber, & Marfatia 2006 (32) 


3.5 


154 


Yuksel & Kistler 2007 133] 


1.7 


221 


Mixed-Iron-Composition: 






Ayq etal. 2005 H 


0.01 


6.6 


Stanev 2008 (351 


0.0002 


1.5 



Kotera, AUard, & Olinto 2010 | 29 | upper 0.08 11.3 

Kotera, AUard, & Olinto 2010 | 29 | lower 0.005 4.1 

Models constrained by Fermi cascade bound: 

Ahlers et al. 2010 (36) 0.09 20.7 

Waxman-Bahcall (WB) fluxes: 

WB 1999, evolved sources (37) 1.5 76 

WB 1999, standard ET) 0.5 27 

In Table |II| we give expected neutrino event totals from a 
wide range of currently allowed cosmogenic neutrino models 
for ARA in three years of operation, compared to recent pub- 
lished expectations for the best current limits to date, from the 
ANITA-II flight |3 1. It is evident that ARA-37 will extend in 
sensitivity above ANITA-2's sensitivity by factors of two or- 
ders of magnitude or more. For strong-source-evolution and 
baseline models, ARA-37 detects between of order 50 to over 
200 events in three years of operation, enough to establish the 
basic characteristics of the energy spectrum and source arrival 
directions. 

There are also recent cosmogenic neutrino flux estimates 
which compute neutrino fluxes subject to constraints from the 
Fermi diffuse gamma-ray background | 36 |, and which include 
a heavier nuclear composition (e.g., an admixture of iron) for 
the UHECRs EHOESl. Over a 3-year timescale afl of these 
models are detectable, but in some cases only marginally, and 
up to five years will be necessary to establish the flux. Over 
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the planned instrument life of a decade or more, ARA-37 will 
thus be able to not only establish the flux levels for all of even 
the most conservative models, but to begin measurements of 
their energy spectral dependence as well. 



C. Resolution 

Although not directly important for detection of neutrinos, 
the resolution of both the distance and angles to the neutrino 
interaction vertex, as well as the ability to reconstruct coarse 
neutrino incident directions on the sky, are important char- 
acteristics of our detector, and we have studied them in de- 
tail. This is especially important for our current realization of 
ARA-37, since the wider spacing will lead to very few multi- 
station coincident events, and thus each station must function 
as a stand-alone neutrino detector in both shower energy esti- 
mation and neutrino direction angular resolution. 

To make these measurements, we have 16 antennas per sta- 
tion, and thus 16 waveform amplitudes and phases, as well as 
the frequency spectral components of the coherently-summed 
waveform which can be estimated to good precision once the 
arrival direction is fitted. From the Vpol and Hpol data we 
also fit the plane of polarization, and with precise timing we 
can measure the radius of curvature of the arriving wavefront. 

Our measurement of the distance to the neutrino vertex is 
accomplished by the estimates of the wavefront curvature. 
This may be thought of as measuring the residuals when fit- 
ting the arrival times to a plane wave. For the angular mea- 
surements, the antenna array is analyzed as a correlation inter- 
ferometer, and precise timing differences between the arrival 
times of the Askaryan radio impulse are determined for all of 
the N{N - l)/2 pairs of N antennas. 

Complementing the precise timing measurements, we can 
also operate our cluster array as a radio intensity gradiome- 
ter and polarimeter. The gradiometric function comes through 
amplitude calibration of the received impulse, and the polari- 
metric information comes from ratios of the calibrated ampli- 
tudes of the Vpol and Hpol antennas. 

All of these estimates are done in offline reconstruction rou- 
tines. They are not necessary for the triggering of the array to 
record potential neutrino events, but they do make maximal 
use of the recorded information in the waveforms and arrival 
times of the events. 



1. Vertex Resolution 

The critical parameter for vertex location is the intra-cluster 
timing precision. For this we have used actual measurements 
made with ANITA data, to which our collaboration has access. 
The ANITA payload, which uses waveform digitizers that are 
comparable to our planned digitizers, has demonstrated tim- 
ing resolution as good as 30 ps rms for waveforms registered 
at the 4a-level detection threshold of ANITA. These timing 
precisions come about from extensive in-flight calibration us- 
ing ground-based impulse generators, and have proven robust 



in the ANITA analysis fS^. For our simulations we have de- 
rated these values by a factor of 3.3 to account for our more 
limited radio bandwidth, the slower sampling rate we expect 
to use, and for possibly unknown systematics in our calibra- 
tion. 



Fig. 28 left,middle) shows the results of these simulations 
for both the range and pointing resolution to the vertex. The 
latter values are important for determining whether an event 
originates under the ice or from above the ice, and the former 
values, combined with our knowledge of the ice attenuation, 
will bear directly on our ability to perform calorimetry on the 
neutrino shower. 



2. Incident neutrino direction resolution 

Estimating the neutrino vertex location in three dimensions 
does not immediately determine the arrival direction of the 
neutrino itself, and it may appear that the tight constraints 
of the cluster geometry would preclude determination of the 
neutrino direction, since this is normally done by imaging of 
the Cherenkov cone, at least in ring-imaging Cherenkov de- 
tectors. However, the richness of the radio wave information 
content does in fact allow for first order Cherenkov cone de- 
termination, using both the amplitude gradient of the cone, 
and the polarization vector, which lies in the plane containing 
the Poynting vector of the radiation and the shower momen- 
tum vector. Thus once the Poynting vector is determined via 
the vertex reconstruction, the plane of polarization combined 
with the local gradient in the cone amplitude is sufficient to 
constrain the neutrino direction on the sky. 

To study the neutrino direction resolution, we have included 
first-order reconstruction algorithms in our Monte Carlo neu- 
trino simulation. These do not yet perform a maximum like- 
lihood minimization which would be very appropriate for this 
complex problem, but instead they perform a S^-grid search 
over variational parameters once the event has been detected 
in the simulation, which includes appropriate thermal noise 
backgrounds. These simulations are computationally inten- 
sive, and involve full 3D ray-tracing of the radio propagation 



through the ice for each tested grid-point. Fig. 28 right) shows 
the results for about 80 detected events in a simulation run 
at £v = 3 X 10^^ eV. The reconstruction for these events was 
about 80% efficient (that is, 20% of the events failed to recon- 
struct), but the reconstructed zenith angle distribution shown 
has a very acceptable standard error of ~ 6°, the simplicity of 
our reconstruction code notwithstanding. 



3. Energy Resolution 

Once the event geometry has been reconstructed, the 
shower energy is determined via standard parametric equa- 
tions B0H42I. For example, a parameterization due to Lehti- 
nen BOl gives the shower energy as a function of radio fre- 
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fractional range error AR/R reconstructed 9 error, degrees Fitted neutrino direction error A0, degrees 



FIG. 28: Left: Distribution of fractional range errors in single-antenna-cluster reconstruction of neutrino interaction vertex. Left inset: the 
distribution of reconstructed vertex range errors vs. range. Middle: Similar to left plot for reconstructed zenith angle of vertex relative to 
antenna cluster. Middle inset: distribution of reconstructed zenith angle vs. range to vertex. Right: Distribution of polar angle errors for full 
reconstruction of the incoming neutrino direction, using vertex reconstruction, amplitude, and polarization information. Right inset: the 2-D 
distribution of reconstructed directions relative to true neutrino direction. 



quency as: 



\Esh{f;RM = V2kR-^iuijo 2 L /sine ^-(^^)'(cose-i/^)V2 

(6) 

where R is the shower distance, Q and L are standard shower 
parameters relating to the excess charge and shower length, 
is the observer's polar angle with respect to the shower (and 
neutrino) direction, k = 27z/X is the wavenumber, n the index 
of refraction of the ice; for typical dielectrics jn= I, and = 
471 X 10~^ is the permeability of free space. 

Using standard variational analysis the fractional shower 
energy (Esh) variation with range error AR is \AEsh/Esh\R = 
e~^^{AR/R) for small ice attenuation coefficient a = L~^. 
The variation with uncertainty in a is \AEsh/Esh\a = 
Re~^^{Aa) which is generally negligible since Aa is of order 
10~^ at most. The corresponding variation with neutrino an- 
gular error A0 is \AEsh/Esh\Q — cos0cA0, where we evaluate 
the variation near the Cherenkov angle 0c. Finally an esti- 
mate of the neutrino energy from shower energy must account 
for the large variance in the Bjorken-j-distribution, defined 
by Ey y~^Esh for charged-current V^,Vx and all neutral- 
current events; electron neutrino events will be less affected 
by this, so our estimate is conservative. At neutrino energies 
of 10^^~^^ eV, numeric evaluations give (y) 0.22, Ay/y 1, 
thus \AEy/Ey\y = Ay/{y) I. Assuming these errors are un- 
corrected, and using AR/R ^ 0.02 with a mean R ^ I km, 
and cosdc 50 = 0.06, the root- sum- squared error is domi- 
nated by the Bjorken-}? uncertainty, giving \AEy /Ey\total ^ 1 
for £'v = 3 X 10^^ eV. This resolution will also be compara- 
ble for lower neutrino energies in the GZK neutrino spectral 
range. The j-dominated uncertainty is generic for UHE neu- 
trino experiments, but this energy resolution is wholly ade- 
quate for the first-order science goals of the ARA instrument. 



D. Comparison to Existing Instruments 

Fig. [29| provides a comprehensive graphic summary of the 
comparison of our estimated ARA sensitivity to estimates for 
several operating experiments, along with 2006 limits from 
the ARA forerunner experiment RICE |2|. We have already 
noted the comparison of ARA to the published ANITA limits; 
here we use projections for ANITA' s reach after three flights, 
along with similar projections for IceCube and the Auger Ob- 
servatory. GZK neutrino models are also included from a wide 
range of estimates |[27l [28l [S^ [HI [SI EZl |43l, including the 
pure-Iron UHECR composition model noted above. 

ARA improves over any other current instrument by an or- 
der of magnitude within 3 years of operation, filling in an im- 
portant gap in sensitivity in the heart of the cosmogenic neu- 
trino spectral energy region. IceCube has excellent sensitivity 
to lower energies, up to the 10 PeV level, and ANITA has un- 
matched sensitivity at the higher energies, above 10 Eev. The 
Auger Observatory, while probing a similar energy range as 
ARA, does not have as high a neutrino sensitivity as it is pri- 
marily a UHECR instrument. ARA will complement these 
other instruments by making high sensitivity observations in 
the 0.1-10 EeV energy range, matching the peak of the ex- 
pected cosmogenic neutrino fluxes. 



V. CONCLUSIONS 

We have described the design and initial performance of 
a new ultra-high energy neutrino detector at the South Pole, 
the 16-antenna, self-triggering ARA-testbed, which is a high- 
fidelity prototype for future ARA detector stations. Our initial 
operation extending well into the the extreme thermal environ- 
ment of the austral winter indicates that radio-frequency inter- 
ference is infrequent and has only a slight impact on operation 
for our testbed detector, which is closest of any future ARA 
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FIG. 29: Compilation of sensitivity estimates from existing instru- 
ments, published limits, and a range of GZK neutrino models, along 
with the expected 3 year ARA sensitivity. 



stations to the primary sources of interference at the South 
Pole station. Other than brief periods of sporadic interference, 
the baseline radio noise levels are dominated by the pure ther- 
mal noise floor of the ambient ice, and the thermal noise does 
not appear to be correlated to wind velocity. We have demon- 
strate the ability to maintain impulse trigger sensitivity at a 
level close to the thermal noise. We have demonstrated RF 
impulse propagation of more than 3 km slant range through 
the South Pole ice without significant loss of signal coherence. 
We have demonstrated inter-antenna pulse timing precision of 
order 100 ps, implying angular resolutions which are more 
than adequate for neutrino vertex reconstruction. We have 
presented simulations using characteristics projected from our 
measurements which give high confidence that our completed 
phase-I array, ARA-37, will achieve its goal of a robust detec- 
tion of cosmogenic neutrinos, and will lay a clear foundation 
for an observatory-class instrument. 
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Appendix A: ARA Autonomous Renewable Power Stations 
(AARPS) 

As ARA moves farther from the station, the transition from 
station power to autonomous power sources will become in- 
creasingly important. The planned ARA footprint calls for 
three ARA stations to be powered from a single node, requir- 
ing about 300W from that node. 

A variety of power sources were reviewed during 2010 in- 
cluding photovoltaic (PV) arrays, wind turbines, diesel gener- 
ators, fuel cells, and Stirling engine generators. The first three 
remain in consideration with the renewable sources, PV and 
wind, to be attempted first. PV is well-known to be rugged and 
efficacious on the Antarctic plateau, so the 2010-2011 season 
was instead used to study wind power options for the Antarc- 
tic plateau. 

The objectives in the first season of deployment were to: 

• Test wind turbine candidates for survivability and for 
power production. While several larger and several 
smaller turbines have been deployed on the plateau, 
there have been no known studies in the intermediate 
power range that ARA requires. 

• Measure wind speed as a function of height from the 
surface. This profile, determined by the turbulence pro- 
duced by the surface upwind of the turbine, will dictate 
the required tower height. 

• Understand tower erection techniques suitable for the 
South Pole environment. 

• Fully instrument the system with sensors for temper- 
ature at various locations, power from turbines, wind 
speeds, etc. connected to a custom System Health Mon- 
itor (SHM). Those data are transmitted to the ICL. 

• Test a Power Instrument Box (PIB). The PIBs deployed 
are thermally insulated but do not have rf shielding; 
shielding is planned for future PIBs. 

• Test for electrical noise in the 200MHz to 800MHz 
band which could result from the turbine itself, the tri- 
boelectric effect on the tower, or emission from instru- 
mentation. 

• Test low-power Zigbee communications. 

The ARA effort to develop autonomous power stations at Pole 
began in the 2010-201 1 season with the deployment of 3 small 
wind turbines, each with a rated power output of slightly over 
a kw: 



At Site 1, a Raum Ikw turbine (Fig. 30) was installed on a 
50-ft lattice tower. It was raised on a tilt-up tower with a Pisten 
Bully and a fixed 10 gin pole. Four guy wires are anchored by 
timbers placed as dead men in trenches in the ice. A single 
anemometer is mounted just below the turbine. 



At Site 2, a Hummer kw turbine (Fig. |3T]) was deployed on 
a 50 monopole. The greater weight of the turbine and tower 
required the use of a Caterpillar 853 to raise the turbine. Three 
sets of four guy wires are anchored by dead men. A single 



22 





FIG. 30: Raum model turbine, installed at South Pole in January, 
2011. 



FIG. 31: Hummer model turbine, installed at South Pole in January, 
2011. 



1. Observations 



anemometer is mounted just below the turbine. A Yagi an- 
tenna is mounted below the anemometer and pointed toward 
the ICL. 



At Site 3, a Bergey XL.l turbine (Fig. 32) was raised on 
a 60 monopole and integral 20 gin pole using a Pisten Bully. 
Three sets of four guys were installed and anchored by dead 
men. In this location only, Kevlar guys are used rather than 
steel cables to reduce static build-up. A Yagi antenna was 
mounted on this tower as well. 

Each turbine was outfitted nearly identically, including the 
following common characteristics, i) A brake switch and 
high-power dump resistors are mounted on each tower, ii) An 



SHM (Fig. [33]) is located inside each PIB to collect metrol- 
ogy and power data and store and/or transmit the data. The 
power is routed to dump resistors inside and outside the box 
(Fig. [34| ), controlled by the SHM to maintain a warmed envi- 
ronment inside the PIB. iii) The SHM provides three destina- 
tions for the data: a 2GB SD card for local storage, a Zigbee 
rf link (operated only for short-term testing until approval is 
obtained), and a hard surface link to the ICL. A Power Board 
containing a three-phase rectifier and measurement of voltage, 
current and rotational rate. 



a. Power Output 

Representative data from one turbine are plotted as a func- 

Very roughly we see the ex- 



tion of wind speed in Fig. 35 



pected relationship of the cube of the wind speed. 



b. Wind Speed as a function of height 

The data collected at three heights on the Bergey tower are 
It shows the height relationship to be 



shown in Figure 36 



relatively mild, but nevertheless illustrates the benefit of taller 
towers, especially at lower wind speed. 



c. Modeling Turbine Performance 

To project the future duty cycle a photovoltaic + wind-based 
system will offer, a detailed turbine simulation has been de- 
veloped. This calculation incorporates turbine- specific man- 
ufacturers' data, in addition to our own 2011 measurements 
of output power vs. wind speed (Fig. 



35 ) and wind velocity 



profile dependence as a function of above-surface elevation 
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FIG. 32: Bergey model turbine, installed at South Pole in January, 
2011. 



(Fig. 36 ). For the particularly high- wind conditions at South 
Pole during the last week of January, 201 1, we find agreement 
between our predicted vs. actual performance at the 5-10% 
level, as shown in Figure [37] Without the aid of batteries, 



based on our model, and folding in measured wind speed data 
from the South Pole over the last half-decade, we expect over- 
all power needs to be satisfied by a photovoltaic+wind-turbine 
system at least 92% of the time. Buffering power with 48-hour 
batteries should improve the duty cycle to >97%. 
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FIG. 33: Schematic of System Health Monitor (SHM) installed at 
South Pole in January, 201 1. 




FIG. 34: System Health Monitor system installed at South Pole in 
January, 2011. 



Power Output vs. Wind Velocity (Hummer Turbine) | 




10 11 12 
Wind Speed (m/s) 



FIG. 35: Plot of turbine output power vs. velocity, illustrating ex- 
pected power law relationship. A fit gives P(velocity)oc v^-^, roughly 
consistent with the expected cubic dependence. 
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I Wind Velocity Ratios as f(Height)~| 

.2 '•'^p 




16.5 m Height Wind Speed (m/s) 



FIG. 36: Ratios of wind speeds at various height combinations, il- 
lustrating the increase in velocity with elevation above snow surface. 



and metrology sensors is operating reliably. Two of the three 
custom power sensor boards continue to operate reliably; the 
one failure may be the Power Board or the cabling and will 
be investigated during the 2011-12 austral summer. The SD 
card operation was verified by retrieving the card after several 
weeks and downloading to a laptop. A simple 4 whip antenna 
was able to receive low-power Zigbee data at the ICL from all 
three sites; a mesh network is planned for next season. Data 
are being received by surface cable using the RS-485 protocol 
from two sites. Received data are being shipped to the north at 
regular intervals. The SHM that was deployed consumes over 
2 watts, but the design cycle and firmware development now 
underway are bringing that value down substantially. A cus- 
tom acquisition and transmission system for the RM Young 
anemometers also operates at very low power. The multiple 
anemometers, temperature sensors and planned vibration sen- 
sor can be daisy-chained to simplify connection. 

e. Power Instrumentation Box 




These were buried in the ice to provide wind protection. 
Although we have some modest temperature control, we do 
not yet have sufficient information to determine how much 
energy is required to maintain the boxes at OC. Structurally the 
PIBs worked well with the proviso that the cable penetrations 
need modification to facilitate access. 



FIG. 37: Comparison of modeled vs. measured wind turbine perfor- 
mance; data taken from Jan. 26-Feb. 1. 



/ Electrical noise from AARP sites 



d. Data Collection 

The SHMs have operated reliably, collecting data from the 
sensors and sending data to the ICL. The array of temperature 



As mentioned elsewhere in this document, attempts to see 
noise from the turbine sites with the test bed have been nega- 
tive. 
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